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Abstract 

We studied the Bouchaud-Mezard (BM) model, which was introduced to explain Pareto's law 
in a real economy, on a random network. Using "adiabatic and independent" assumptions, we 
analytically obtained the stationary probability distribution function of wealth. The results shows 
that wealth-condensation, indicated by the divergence of the variance of wealth, occurs at a larger 
J than that obtained by the mean-field theory, where J represents the strength of interaction 
between agents. We compared our results with numerical simulation results and found that they 
were in good agreement. 
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I. INTRODUCTION 



Researchers in the field of complex networks agree that a change in the network topology 
induces a critical change in dynamics. Pastor-Satorras and Vespignani first showed the 
absence of an epidenhc threshold in a S ea,e- fre e netwo^, foHowin g which m anv res ea r ehe re 
have focused on the dynamics on complex networks, such as synchronization jj, [J, pattern 
formation^], and other phenomena. 

, n 

In this study, we focus on the Bouchaud-Mezard(BM) model on a complex network [5|]. 
It is known that the wealth distribution in a real economy exhibits a power-law behavior, 
called Pareto's law|6i]. With a view to this power law, Bouchaud and Mezard proposed a 
model, given by the following Stratonovich stochastic differential equation. 



where Xi, N, J, and a 2 represent the wealth of the z-th agent, number of agents, coupling 
between agents, and variation of noise, respectively. In this model, the evolution of the 
wealth is determined by two processes: exchange of wealth and a random multiplicative 
process, respectively described by the first and second term in the right-hand side of Eq.(JT]). 
Bouchaud and Mezard analyzed this model using the mean-field theory and calculated the 
probability density function(PDF) of wealth. They showed that the stationary distribution 
of normalized wealth Xij (x), where (• ■ ■) represents the average over all agents, exhibits the 
power-law behavior. They also found that "wealth-condensation," which is indicated by 
the divergence of the variance of Xi/(x), occurs at J < J c = a 2 . The divergence of the 
variance implies that wealth condenses to a few rich agents. On the other hand, if J > J c , 
the variance remains finite, and the wealth of many agents is close to the average. 

In the original BM model, all agents are coupled with each other. However, in a real 
economy, agents can exchange their wealth with a limited number of agents. Therefore, it 
is natural to extend the BM model on a complex network in which the number of neighbors 
is limited. Some studies have already dealt with this subject. In their original study on the 
BM-modelfsj], Bouchaud and Mezard carried out numerical simulations on a regular random 
graph to estimate the exponent of the power-law behavior. They reported that the exponent 
obtained from numerical simulation becomes smaller than that obtained from mean-field 
theory. Some studies have also reported on the simulation of this model on a Barabasi- 
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Albert (BA) network and a Watts-Strogatz(WS) networkj?, 8]. In these studies, the authors 
discussed numerical simulation results, however, none of them proposed a quantitative theory 
that could explain these results. 

This study aims to develop a quantitative theory of the BM model on a random net- 
work. The key assumption of our theory is "adiabatic and independent" assumptions on 
the stationary distribution function, which is explained in a later section. By using these 
assumptions and the central limit theorem, we analytically derive the equations that deter- 
mine the stationary distribution function in the non-wealth-condensate phase. We compared 
our analytic results with those of the numerical simulation, and we found that our theory 
showed better agreement than did the mean-field theory. 

The remainder of this paper is organized as follows. In the next section, we define the 
model we investigate in this paper. Then, we describe our theory and its results, first 
for the case of a regular random network and then for a random network with arbitrary 
degree distribution. In Sec. IIV[ we describe a comparison of our analysis and the numerical 
simulation. In the last section, we summarize our results and discuss the problem to be 
solved. 



II. MODEL 



The original BM model is expressed using the Stratonovich stochastic differential equa- 
tion; however, we use the equivalent Ito stochastic equation for mathematical convenience. 
We consider the BM model on a complex network described by the following Ito stochastic 
differential equations 

N 

)]dt + V2axidW i: (2) 

i=i 

where Xi, J, N, and a 2 are the same as Eq.(TjQ), and represents the adjacent matrix. On 
the network model, we consider a random network in which the degree distribution is given 
by Q(k). 
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III. THEORY 



In this section, we consider the stationary PDF of the normalized wealth of Eq.(j2]). 
Eq.()2]) is invariant under the change of scale x\ = axi for any positive constant a, and we 
can assume (x) = 1 without loss of generality. We derive the analytic form of the stationary 
PDF of Xij (x), the normalized wealth at node i. First, we explain our method on the regular 
random graph, in which all nodes have the same degree k. Then, we extend the analysis to 
the general random network model, whose degree distribution is given by Q(k). 



A. Case of a regular random network 



This subsection focuses on the analysis of the system when each node has the same degree 
k, in other words, Q{k) is a delta function. We assume that pi(x,t), the PDF of wealth at 
node i, is independent of i, Pi(x,t) = p(x,t). 

First, we review the mean-field treatment of the BM model. We consider the distribution 
of wealth at node i. By using the mean-field theory, we approximate | ^ . a^Xj = (x) = 1 
in Eq. ([2]). Under this approximation, p(x,t) satisfies the following Fokker-Planck equation, 



dp 
dt 



and we find p^ q (x), the stationary PDF obtained from this equation, as 



Peq F ( x ) = Cexp(— a/x)x 



-2-a 



(3) 



(4) 



where a = Jk/a 2 and C = a 1+Q /r(l + a). In this case, wealth-condensation, defined as the 
divergence of (x 2 ), occurs at Jk < a 2 . 

To proceed beyond the mean-field approximation, we make what we call "adiabatic and 
independent" assumptions. We define the "local" field Xi = }:J2j a ij x j- Because is not 
only if nodes i and j are connected, Xi represents the local average of the wealth around 
node i. 

If Xi is fixed, the PDF of x at node i is given by solving the " local" Fokker-Planck equation 



dp 9 ri ,,. , 2 n 2 d 



and the conditional PDF of x under the local field x is given by 

p eq (x\x) = C(x) exp(— ax/x)x~ 2 ~ a , 



(5) 



(6) 



where C(x) = (ax) 1+a /T(l + a). 

Here we make the "adiabatic and independent" assumptions. First, we assume the static 
PDF p eq (x) can be approximated by 

Peq(x) = J dxP(x)pi(x\x), (7) 

where P(x) is the PDF of x. If x, changes much slower than x, this condition is satisfied, 
and therefore we call it the "adiabatic" assumption. Under this assumption, the problem to 
calculate p eg (x) is reduced to the one to calculate P(x). 

The second assumption is needed to calculate P(x). We assume that the random variables 
Xj, where j runs in the neighborhood of node i, are independent. This assumption enables us 
to use the central limit theorem to obtain P(x). We consider the case in which the variance 
of x is finite, and that the average and variance of x are 1 and s 2 , respectively. Using the 
central limit theorem, we can approximate 

/ k(x-l) 2 \ 

P(x) ~ exp - 1 - - } . (8) 



Inserting Eq. (jSJ) into Eq. (j^j), we obtain 



Peq(x) ~ / dx j]L exp (~ k ^ X 2 ^ ) peq(x\x). (0) 

v 2ns V £s 



o 



The final step is to check that (x) = 1 and to calculate s 2 . Using J °° dxp eg (x\x)x = x 
and J °° dxp eq {x\x)x 2 = -^tix 2 for a > 1, we obtain 

/poo poo 
xp eq (x) = / dxP(x)x ~ / xP(x)dx = 1, (10) 
JO J -oo 



and 



.2 ... _ / J™ a z.2 ■ a ^ 



x 2 p eq = dx -i 2 P(i) 1 + - . (11) 

1 a — 1 a — IV k 



Here we change the lower limit of integration from to — oo, assuming k/s 2 to be large. 
From Eq. dlOp . we conclude that (x) — 1 , and that there is no inconsistency. From Eqs. ( FlOl) 
and (TIT]) , the variation of x is given by -^r(^~ + !)■ Therefore we find the self-consistency 
condition for s 2 as 

[(A; - l)a - A;]s 2 = fc. (12) 
Eqs.([6]), OH]), and ( fl2l) are the set of equations that determine p eq (x). 
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We should make several comments on these results. First, this theory gives a distribution 
that does not obey the power- law behavior. Eq.Q shows that the PDF p eq (x) is written 
with the integral of P(x)p eq (x\x) over x. Because P(x) is Gaussian, p eq (x) does not exhibit 
the simple power-law behavior. 

The next important suggestion of this theory concerns the wealth-condensation transition. 
In our theory, s 2 diverges if (k — l)a — k = 0, implying that the wealth-condensation occurs 
at this point. Using a = Jk/a 2 , this leads to the conclusion that wealth condensation 
occurs at J c = a 2 / (k — 1). On the other hand, the mean-field theory gives the divergence of 
the variance at J c = a 2 /k. Therefore the difference between our theory and the mean-field 
theory can be tested by estimating J c from the simulation. 

Finally, we note that our theory can be applied only for the non-wealth-condensation 
phase. We need the central limit theorem to obtain P(x), which is only applicable for the 
case in which the variance of x is finite. We discuss this point in greater detail in the final 
section. 



B. Case of a general random network 

In this subsection, we extend the developed method for a regular random graph to a 
general random network. 

As in the case of a regular random network, we start from the mean-field theory. In this 
theory, we define pk{x,t) as the PDF of x on the node whose degree is k. Because this 
network is heterogeneous, the mean of x may depend on k. Therefore, we need to define 
Xk = J dxxpk{x,t), the average of x on nodes with degree k, to perform the mean- field 
calculation. The mean-field Fokker-Planck equation is constructed in the same manner as 
in the case of the SIS model or Kuramoto oscillator , and we obtain 



dp k (x,t) 
dt 



d_ 

dx 



kk Q(k) ,_ \ i ±\ 
'x k > - x)p k {x,t) 



(k) 



— a xpk{x, t) I + a 



d_ 

dx 



x-^(xp k (x,t)) 



(13) 



We can simplify this equation by introducing the "weighted" average of x as x = 



kQ(k) 
k (k) 



Xki 
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which leads to 

(14) 



dpk d 2 2 <9 

- -— [Jfc(x - xjpfc - a xp k \ + o 



x^(x Pk ) 



dt dx ' dx 

where we abbreviated pk(x,t) as pk- As in the case of a regular random graph, we obtain 
Peqjt( x )i the stationary PDF by the mean-field theory, as 

pf q F k (x) = C k (x)exp(-a k x/x)x- 2 - ak , (15) 

where a k = Jk/a 2 and C k (x) = (akx) 1+ak /T(l + oik)- Finally, x is determined to satisfy the 
condition (x) = 1. From f dxxp^Xix) = x, we obtain x = 1. 

Now we follow the same procedure as that described in the previous subsection to calculate 
the PDF more accurately. We consider the PDF of the wealth on node i, whose degree is 
k. If the average of sums of x in the neighborhood of node i is given by x, we find the 
conditional PDF of x on node i as 

pk{x\x) = Ck(x)exp(—akx/x)x~ 2 ~ ak , (16) 

Using the adiabatic assumption explained in the previous subsection, we assume 

Pk{%) = J dxP k (x)p k (x\x). (17) 

Next, we approximate Pk{x) by Gaussian using the independent assumption. One dif- 
ference between the regular random graph and the random graph with arbitrary degree 
distribution lies in that we need to assume Lindeberg's condition in this case. Suppose that 
there are independent variables yi,---y m , whose mean and variance are pi, P2, • • • Pm an d 
sf, s|, • • • respectively. Then, the PDF of — (yi + y 2 + ■ ■ ■ y m ) converges to the Gaussian, 
whose mean and variance are £*2 + ' ■ ■ + Pm) and —s(cr 2 + o\ + • • • <^ 2 n )-, for m — > oo, if 



Lindeberg's condition is satisfied 10j] . Although m is finite in our case, we can approximate 



where S% is a parameter that must be determined by the self-consistency condition. 

The final step is to obtain self-consistent equations for the mean and variance of x k - 
Using J dxxpk(x) = 1 and f dxx 2 pk(x) = + S%) obtained from Eqs. flTB"]) . (IT7|) and 

( JT8l) . we respectively obtain Hk and s|, the average and variance of x on nodes with degree 
k, as pk — 1 an d 

4 = £*S+1. d9) 

ak - 1 



A node with degree k is connected to a node whose degree is k' with probability kk ^ k - , 
and we obtain 

1 x - kk'Q(k') 2 u 

Sk = v\?^r Sk ' = v (20) 

where u represents the weighted average of the variance, 

u = J2kQ(k)s 2 k /(k). (21) 

k 

From Eqs. (JT9|) and (l2"Uj) . we obtain 



n a^u/k + 1 

CK fc - 1 

Inserting this equation into Eq. (l2"Ti) . we obtain 

Q(k) ctkU + k 

k 

which leads to our final result, 



(22) 



^ (k) a k - 1 ' 1 ; 



\- Q{k)a k \ \- Q(fc)fc 

Eqs.f|T6l). ( |T7I) . ( fT8i) . ( 1201) . and ( 1241) are the set of equations that determine the stationary 
PDF. 

The most important difference between the mean-field theory and ours lies in that there 
is no effect of Q(k) in the mean-field theory. In the mean-field theory, the PDF on the 
node with degree k given by Eq. (fl5|) does not depend on Q{k), which causes a curious 
behavior. Suppose that we decrease J from a large value to 0. At J = <r 2 /k, the node with 
degree k goes into the wealth-condensate phase, whereas nodes with a larger degree do not 
condensate. Therefore, the mean-field theory predicts the coexistence of condensated and 
non-condensated nodes. In this theory, low-degree nodes condensate at large J, whereas 
high-degree nodes do not condensate until J becomes sufficiently small. 

On the other hand, all PDFs with different degrees are connected through u in our theory. 
In this theory, u diverges when ^ k m^^zu = 1> which implies the divergence of all Sk from 
Eq.( )22|) . In our theory, wealth-condensation occurs on all nodes simultaneously. 

Finally, we comment on the behavior of the BM model on a scale-free network. Re- 
searchers in the field of complex networks may consider that a singular behavior occurs 
in a scale-free network upon the wealth-condensation transition, such as the divergence or 
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disappearance of the transition point J c . Unfortunately,this is not the case. The main dif- 
ference between our case and other models that show singularity, such as the SIS-model or 
Kuramoto transition, lies in that we impose the self-consistency condition on the variance 
of x, and not on its average. As shown in Eq. (EOT) . we divide the weighted average of the 
variance s 2 . by k to calculate Sf.. This eliminates the singular behavior that we often observe 
in the dynamics on a scale-free network. 

IV. SIMULATION 

In this section, we test the analytic results obtained in the previous section by comparing 
them with numerical simulations. 

Because we have presented two analytic results, one for a regular random network and the 
other for a heterogeneous random network, we carried out the simulations for both networks. 
For the former, all nodes had the same degree, whereas for the latter, half of the nodes had 
degree 10 and the other half had degree 20. 

We use these two models because it is easy to calculate the PDF in these models. However, 
the readers might find these models too artificial. To test our theory in a more realistic 
model, we show the result of simulations on B A- network [11]. Here we note that the direct 
calculation of wealth distribution is rather hard because we must calculate integral in Eq. (fT7|) 
for many values of k. In this paper, we investigate s 2 ., the variation of x at a node with 
degree k to test our theory instead, because it is much easier to compute. 

In the following simulation, we carried out numerical integration by the Euler-Maruyama 
algorithm to obtain the distribution of normalized wealth x/(x). 

A. Case of a regular random network 

In this subsection, we show the simulation result for a regular random network and 
compare it with our theory. 

First we calculate the PDF of the normalized wealth from 10 simulation trials on networks 
that include 5000 nodes. In FigfTJ we show the obtained PDF for the network with degree 
k = 10, a 2 = 1 and coupling J = 0.3, 0.4, and 0.5. For all J, the distribution coincides 
well with our theory, indicated by the solid line. The PDF given by the mean-field theory, 
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N = 5000, k = 10, J = 0.3 
PDF 



N = 5000, k = 10, J = 0.4 
PDF 




X/<X> 1.x 10-4 



N = 5000, k = 10, J = 0.5 
PDF 




x/<x> 




1.x 10 



x/<x> 



FIG. 1: Log-log plot of PDF of normalized wealth obtained by simulation when k = 10, a 1 = 1, 
and J = 0.3, 0.4, and 0.5. The solid and dashed lines indicate the result obtained by our theory 
and that obtained by the mean-field theoryrespectively. 
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FIG. 2: Temporal and spatial behavior of x. Left: time series of x and x at one node. Right: 
scatter plot of x and x. Both plots are obtained from the simulation at N = 5000, k = 10 and 
J = 0.3. 

indicated by the dashed line, strongly underestimates the probability density at small x/(x). 
On the other hand, our theory gives a slightly larger PDF than the simulation, although, it 
shows better agreement. 

Here we investigate the temporal and spatial correlation of x to test the "adiabatic and 
independent" assumptions used to derive Eqs.flZj) and j9j). In the right-hand side of Figj2j 
we plot the typical trajectory of x{t) and x{t). Because the model we study is described 
by stochastic differential equations, it is difficult to show that the change in Xi is "slower" 
than that in Xi, however, it seems that xi changes more quickly than X{. Concerning the 
"independent" assumption, we show the scatter plot between x and x in the left-hand side 
of FigJ21 Although there exists a tendency for x to increases as x increases, it is not strong. 
The correlation between x and x calculated by the numerical simulation is 0.238. Therefore 
we conclude that the "adiabatic" "independent" assumptions are fairly good. 
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N = 5000, k = 10, J = 0.5 N = 5000, k = 20, J = 0.5 

PDF PDF 




FIG. 3: Log-log plot of PDF obtained by simulation when J = 0.5, a 2 = 1 and k = 4, 10, and 20. 
The solid and dashed lines indicate the result obtained by our theory and that obtained by the 
mean-field theory, respectively . 

In Fig. [31 we show the PDF at J = 0.5, a 2 = 1 and k — 4, 10, and 20. It is clear that the 
discrepancy between our theory and the numerical simulation increases as k decreases. This 
is not surprising, because the error caused by the application of the central limit theorem 
increases at small k. However, our theory appears to be better than the mean-field theory, 
even in the worst case. For example, when k = 4, the probability density obtained by the 
mean-field theory is smaller than 10~ 5 at x — 0.1, whereas the numerical simulation and 
our theory shows that it is O(10 _1 ). On the other hand, the difference from the result 
of numerical simulation at large x is indistinguishable between our theory and mean-field 
theory. 

Finally, we check the wealth-condensation transitions. Though our theory does not ex- 
hibit the exact power-law behavior, it suggests that variance of normalized wealth diverges 
at J c = o 2 j(k — 1). Therefore, we can expect power law behavior p(x) oc x~ 3 for large x 
at this point. On the other hand, the mean-field approximation gives J c = cr 2 /k, and we 
can test our theory by estimating J c by checking the exponential tail. In Fig|4]we plot the 
exponent 7, p(x) oc x -7 at x 3> 1, estimated from the fitting of the PDF obtained through 
100 simulation trials on a regular random graph iV = 1000, k = 4. It is clear that wealth- 
condensation occurs at J ~ 0.35, which is closer to the prediction of our theory J c = 0.33, 
than that of the mean-field theory , J c = 0.25. 

B. Case of a heterogeneous network 

In this section, we present the simulation result for a heterogeneous network. As explained 
in Sec JIIIl the PDF p eq ,k(x) does not depend on Q{k) in the mean-field approximation. In 
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FIG. 4: Exponent obtained by simulation for k = 4, a 2 = 1. The dotted line indicates the 
wealth-condensation transition. 

this approximation, the PDF of the wealth for nodes whose degree is k equals that for the 
regular random graph. On the other hand, our theory predicts that the PDF changes if 
Q(k) differs. To test our theory for heterogeneous network, we make the network that half 
of the node have degree ki = 10, and other half have degree k 2 = 20. 

The PDF of the normalized wealth obtained by the simulation is shown in Figj5] for 
J = 0.3 and a 2 = 1. In this figure, we also plot the PDF obtained by our theory for a 
heterogeneous network and a regular random graph with degree k\ or k 2 and by the mean- 
field theory and that by our theory for regular random graph with degree k\ or k 2 . For the 
PDF on nodes with degree 10, the PDF obtained by our theory, indicated by the solid line, 
is suppressed at small x compared with the case of a regular random network, indicated 
by the dashed line. The PDF obtained by the simulation shows better agreement with our 
theory for a heterogeneous network. For the PDF on node with degree 20, the difference 
among the three theories is small; however, we can conclude that our theory can estimate 
the PDF very well. 

C. Case of B A- network 

To test our theory in a more realistic network model, we show the simulation result on 
the BA-network model. The PDF obtained from simulations on a BA-model with minimum 
degree 4, J = 0.3 and o 2 = 1 is given in the left-hand side of Figj6j As we have noticed, 
it is slightly difficult to show the PDF obtained from our theory. Instead of the PDF, we 
plot the s 2 ,, the variation of x on nodes whose degree is k, obtained by simulation and that 
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FIG. 5: PDF obtained from numerical simulation on a heterogeneous network. The solid, dashed, 
and dotted lines indicate the PDF obtained by our theory for a heterogeneous network and a 
regular random graph, and by the mean-field theory, respectively 




FIG. 6: Result of simulation for BA-network with minimum degree 4 and J = 0.3, a 2 = 1.0. Left: 
PDF. Right: s?, the variation of x on nodes with degree k. Thick line indicates the theoretical 
result obtained from Eq. (|22p . 

obtained by Eq. in the right-hand side of Fig. [5]. The theoretical prediction, indicated 
by the thick line, coincide well with the result of numerical simulation. 



V. SUMMARY 



In this paper, we derived the PDF of the Bouchaud-Mezard model on a random network. 
Using adiabatic and independent assumptions, we derived the equations for the PDF, Eqs. 
(ED, (ED, and (H2D for a regular random network, and Eqs. (DID, (QTfl, (JT8|). (|20p. and (gH) 
for a general random graph. It is difficult to solve these equations analytically; however, 
doing so provides considerable information about the wealth-distribution. In particular, 



13 



we can analytically obtain the wealth-condensation point J c for a regular random graph. 
These analytic results are compared with the numerical simulations in Sec. IIVt and good 
agreement is found with our theory. Below, we discuss the problem to be solved. 

First, we note that the approach we have used in this paper cannot be applied to a wealth- 
condensate phase. We use the central limit theorem to derive the equation for a stationary 
PDF; however, this approach is not applicable when the variance of p(x) diverges. To treat 
the wealth-condensate phase, we need a generalized central limit theorem that is applicable 
when the variance diverges. It is known that the PDF of the sum of independent identical 
random variables converges to a stable distribution function when the variance diverges 10j. 
It would be possible to make the theory for the wealth-condensate phase on a regular random 
network using this theorem; however, further work will be required toward this end. For a 
general random network, the generalized central limit theorem unfortunately remains still 
insufficient. To treat such a network , we require a generalized central limit theorem for 
non-identical independent random variables, which has not been yet established. 

Second, one might be interested in better approximating the the PDF. As shown in Fig. |3j 
the PDF obtained by our method still has a large discrepancy with the numerical simulation, 
especially for nodes with small degree. To reduce these errors, we need to use the improved 
central limit theorem, which can deal with the rate of convergence. The application of 



Chebyshev's rate-of-convergence theorem 12| would lead to a better approximation. 

Finally, we should comment on the "adiabatic and independent" assumptions. We have 
no proof to verify these assumptions; however, the coincidence with the simulations suggests 
that these assumptions work very well. If so, the theory developed in this study will be 
applicable to other dynamical systems. The essential part of our theory is to impose the self- 
consistency condition on the average and variance of the dynamical variables on each node. 
In the case of the BM model, the self-consistency condition for the average is automatically 
satisfied, and we need only to calculate the variance. This procedure is very general, and it 
will be applicable to analyze other dynamical behaviors such as synchronization or diffusion. 
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